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ABSTRACT 

Supermassive binary black holes (BBHs) are unavoidable products of galaxy mergers and are expected to 
exist in the cores of many quasars. Great effort has been made during the past several decades to search 
for BBHs among quasars; however, observational evidence for BBHs remains elusive and ambiguous, which 
is difficult to reconcile with theoretical expectations. In this paper, we show that the distinct optical-to-UV 
spectrum of Mrk 231 can be well interpreted as emission from accretion flows onto a BBH, with a semimajor 
axis of ^ 590 AU and an orbital period of ^ 1.2 year. The flat optical and UV continua are mainly emitted 
from a circumbinary disk and a mini-disk around the secondary black hole (BH), respectively; and the observed 
sharp drop off and flux deficit at A ~ 4000 — 2500A is due to a gap (or hole) opened by the secondary BH 
migrating within the circumbinary disk. If confirmed by future observations, this BBH will provide a unique 
laboratory to study the interplay between BBHs and accretion flows onto them. Our result also demonstrates 
a new method to find sub-parsec scale BBHs by searching for deficits in the optical-to-UV continuum among 
the spectra of quasars. 

Subject headings: accretion, accretion discs - black hole physics - galaxies; active - galaxies: nuclei - galaxies: 
individual (Mrk 231) - quasars: supermassive black holes 


1. INTRODUCTION 

Supermassive binary black holes (BBHs) are natural prod¬ 
ucts of the hierarchical mergers of galaxies in the ACDM cos¬ 
mology and are expected to be abundant (e.g..lBegelman et al.l 
[Tom iMerritt & Milosavlievidl2005h . since many 

galaxies (if not all) are found to hos t a supermassive black 
hole (SMBH) at their centers (e.g.. iMagorrian et"^ 119981: 
iKormendv & Hoir2013h . Evidence has been accumulated for 
SMBH pairs in active galactic nuclei (AGNs) and quasars 
with perturbed galaxy morphologies or other merger features 
(e.g.,jKomossa et al 120031 : [Liu et alJl2010l : IComerford et alJ 
1201 U lEu et al.ll2012 ). These SMBH pairs will unavoidably 


evolve to closely bound BBHs with separations less than 1 pc. 
However, th e evidence for BBHs at the sub-pc scale is still 
elusive ('e.g.. lPopovicll201^ . which raises a challenge to our 
understanding of the BBH merger process and the formation 
and evolution of SMBHs and galaxies. 

A number of BBH candidates in quasars have been pro¬ 
posed according to various spectral or other features, such 


sion (e.g.. IBoroson & Lauerl 

o 

o 

Tsalmantza et all 1201 It 

lEracleous et alJI2012l lJu et al. 

20131 

Liu et al.ll20l4). the ne- 


120151) . and etc.; however, most of those candidates are still 
difficult to be confirmed. Thus, it is of great importance 
to find o ther ways to select and i dentify BBHs in quasars. 
Recently, iGiiltekin & MilleJ (1201 2^ proposed that the contin¬ 
uum emission from a BBH-disk accretion system, with unique 
observable signatures between 2000A and 2/rm because of 
a gap or a hole in the inner part, can be used to diagnose 
BBHs (s ee ISesa na et al.l 1 201 2t iRoedig et al1l2014t lYan et al.l 
12014 but lEarris et al.ll2015 ). This method may be efficient in 


identifying BBHs since many AGNs and quasars have multi¬ 
wavelength observations and broad band spectra. Those pre¬ 
vious investigations only focus on theoretical predictions, and 


the present paper is the first attempt to apply this method to fit 
real observations. 

In this paper, we report a BBH candidate in the core of 
Mrk 231, the nearest quasar with a redshift z = 0.0422, ac¬ 
cording to its unique optical-UV spectrum. In Section |2] we 
summarize the multi-wavelength spectrum of Mrk 231 and its 
distinctive spectral features comparing with normal quasars. 
The spectrum of Mrk 231 at the optical band is similar to the 
quasar composite spectrum; however, it drops dramatically 
at the wavelengths around 3000A and becomes flat again at 
< 2500A. This anomalous contin uum is hard to be exp lained 
by normal extinction/absorption (IVeilleux et al.l |20T^ . We 
propose that the unique optical-to-UV spectrum of Mrk 231 
can be explained by emission from a BBH accretion system, 
with which the drop of the continuum at < 4000A is due to 
a gap or a hole opened by the secondary component of the 
BBH. In Section [3 we introduce a simple (triple-)disk model 
for the accretion onto a BBH system. Using this model, we fit 
the optical-to-UV continuum of Mrk 231 (SectionlDi and con¬ 
strain the orbital configuration of the BBH system and the as¬ 
sociated physical parameters of the accretion process in Sec¬ 
tion |5] Discussions and conclusions are given in Sections |6] 
and|2] 


2. MULTI-BAND OBSERVATIONS OF MRK 231 

Mrk 231 is an ultraluminous infrared galaxy with a bright 
quasar-like nucleus. It is probably at the final stage of a 
merger of two galaxies as suggested b y its disturbed mor¬ 
phology and the associate d tidal features (iLinari & Macchett^ 
ll994TArmus et akl 119941) . The broadband spectrum of the 
Mrk 231 nucleus exhibits some extreme and surprising prop¬ 
erties as follows. 

Eirst, the flux spectrum {F\) drops dramatically by a factor 
of ~ 10 at the near UV band (from wavelength A ~ 4000 A 
to 2500 A), while it is flat at A ~ 1000 A- 2500 A and at A ^ 
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4000 A-10000 A. If this sharp drop off is due to extinction, it 
requires a large dust reddening of Av ~ 7 mag at A ~ 2500 — 
4000 A and a small dust reddening ^ 0.5 mag at < 2500 A 
(IVeilleux et al.llMoh . Thus, an unusual complex geometric 
structure for the extinction material surrounding the di sk must 
be designed (IVeilleux et al.ll201^lLeighlv et al]l2014li . 

Second, Mrk 231 is an extremely powerful Fe II emitter as 
suggested by the spectrum on the blue side of Hq, and on both 
sides of Such a high level of optical Fe line emission 
suggests a significant amount of Fe line emission at the UV 
band; however, it is not visible in the ob served UV spectrum 
(IVeilleux et al.llMT^lLeighlv et al.ll2014ft . 

Third, a number of broad low-ionization absorption line 
systems, such as Na I D, Ca II, Mg II, Mg I, and Fe II, have been 
identified in the optical-to-UV spectrum, which suggests that 
Mrk 231 should be a n Fe low-ionization b road absorption line 
quasar (FeLoBAL) (IVeilleux et al.ll20T3t) . However, the ex¬ 
pected corresponding absorption features at the UV and FUV 
bands are not evident. 

Fourth, the hard X-ray emission of Mrk 231 is extremely 
weak. The intrinsic, absorption corrected. X-ray luminosity 
at 2-lOkeV is L 2 -iokeV ~ 3.8 x lO'^^erg s“^, and the ra¬ 
tio of L 2 -iokeV to the bolometric luminosi ty (inferred from 
the optical lumin osity) is only ^ 0.0003 (ISaez et al.l 120121 
iTeng et akl 120141) . almost two orders of magnitude smaller 
than the typical yal ue 0.03) of a quas ar with a similar bolo¬ 
metric luminosity (iHopkins et al.ll200^ . 

Below we show that the aboye first feature, the anoma¬ 
lous UV continuum, is a distinct prediction of a BBH-disk 
accretion system as shown in Figure [U The last three fea¬ 
tures, commonly interpreted under the context of complicated 
outflows with absorptions, can be also accommodated under 
the framework of the BBH-disk accretion, e.g., the wind fea¬ 
tures are consistent with the Fe absorption features of a typical 
FeLoBAL, and Mrk231’s intrinsic X-ray weakness is also a 
natural consequence of a BBH-disk accretion system with a 
small mass ratio. 

3. OPTICAL-TO-UV CONTINUUM FROM A BINARY BLACK 
HOLE—(TRIPLE-)DISK ACCRETION SYSTEM 

Considering a BBH system resulting from a gas rich 
merger, the BBH is probably surrounded by a circumbinary 
disk, and each of the two SMBHs is associated with a mini¬ 
disk (see Fig. [T]i. In between the circumbinary disk and 
the inner mini-disks, a gap (or hole) is opened by the sec¬ 
ondary SMBH, which is probably the most distinct feature 
of a BBH-disk accretion system, in analogy to a system 
in which a gap or hole is open ed by a planet migrating in 
the pl anetary disk around a star (iLin et al.ll 19961: lOuanz et al.l 
120131) . This type of geometric configurations for the BBH- 
disk accretion systems has be en reyealed by many numer- 


ical simulations and analvsis (lArtvmowicz & Lubowl 

1996 

iDorazio et alJl2flL3b iHavasaki et al.1 

12008b iFscala et all 

2005 


ICuadra et al.1 l200l: iRoedig et al.112^ iFarris et all |2^ P1 


The continuum emission from disk accretion onto a BBH 
may be much more complicated than that from disk accre¬ 
tion onto a single SMBH, since the dynamical interaction be¬ 
tween the BB H and the accretion flow onto it changes the 
disk structure (iGiiltekin & Milled 120121; ISesana et akl 120121; 

* The width of the gap (or hole) is roughly determined by, but could be 
somewhat lai'ger than, the Hill radius Rn. However, set a slightly large gap 
size, e.g., does not affect the results presented in this paper signifi¬ 

cantly. 



Fig. 1.— Schematic diagram for a BBH-disk accretion system. The BBH 
is assumed to be on circular orbits with a semimajor axis of asBR. and the 
masses of the primary and secondary components are M#,p and re¬ 

spectively. The BBH is surrounded by a circumbinary disk, connecting with 
the mini-disk around each component of the BBH by streams. In between 
the circumbinary disk an d the inner mini-disks, a gap or hole is opened 
by the secondary SMBH jArtvmowicz & LubQwlll99q;IDorazio et alJI2Q13l; 
IFarris et '^120141) . The width of the gap (or hole) is roughly determined by 
the Hill radius Rn ~ 0.69Q^/^aBBH]-where 

q is the mass ratio, and the inner boundary of the circumbinary disk can be ap¬ 
proximated as Tin c ~ ^BBH / (1 +?)+-Rh • The outer boundary of the mini¬ 
disk surrounding the secondary SMBH (rout,s) is assumed to be a fraction 
(/r,s) of the mean Roche radius, -Rrl(^) — 0.49aBBHQ^^^/[0.6g^/^ + 
ln(l + <Eggleton| f l983f) , i.e., rout,s = /r,s-RRL(g), considering 

that the mini-disk may not fill the whole Roche lobe (the red dashed cir¬ 
cle). For BBHs with mass ratios roughly in the range from a few percent to 
0.25, the accretion onto the secondary SMBH and consequently its emission 
dominates, compared with that from the mini-disk around the primary BH 
IRoedig et alJ2()l3 ; IFarris et alJ2Q14l) . 

Raflkovl2013l ; lRoedig et al.l2014l ; [Yan et al.l20141; IFarris et alJ 

20151) . Nevertheless, we adopt a simple model to approx¬ 
imate the continuum emission from a BBH-disk accretion 
system as the combination of the emissions from an outer 
circumbinary disk and an inner mini-disk around the sec¬ 
ondary SMBH, each approximated by m ulticolor black body 
radiation in the standard thin di sk model (IShakura & Sunva^ 
ll973l:lNc)vikov & Thorne! 19731) . The emission from the mini¬ 
disk around the primary SMBH is insignificant for a BBH 
system with a small mass ratio (roughly in the range of a few 
percent to 0.25) due to its low accretion r ate as suggested by 
the state of the art numerical simulations (iRoedig et al.ll2012l ; 
IFarris et aP 120141) . thus its emission can be neglected. Our 
analysis suggests that a large q cannot lead to a good fit to the 
observations. 

3.1. Emission from the circumbinary disk 

We choose a standard thin disk to approximate the tempera¬ 
ture profile of the circumbinary disk. The structure and spec¬ 
tral energy distribution (SED) of the circumbinary disk may 
be different from that of a standard thin disk, especially at 
the region close to the inner edge, since the torque raised by 
the central BBH may lead to gas accumulation there. In this 
region, the circumbinary disk is somewhat hotter than the cor¬ 
responding region of a standard thin disk with the sa me accre- 
tion rate and total SMBH mass (M, p + M, s) (e.g.. lRaflkovl 
120131) . We neglect this slight difference for now and will dis¬ 
cuss this in the Appendix. 
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3.2. Emission from the inner mini-disk associated with the 
secondary SMBH 

We assume that the emission from the inner mini-disks is 
dominated by that from the mini-disk around the secondary 
SMBH, and the emission from a mini-disk can also be ap¬ 
proximated by that from a standard thin disk with the same 
extent, accretion rate, and SMBH mass. The a ccretion onto 
each of the two SMBHs ma y be highly variable (iRoedig et al.l 
120121 iHavasaki et al.ll2()()^ . and the temperature structure of 
each of the two mini-disks may be affected by the torque 
from the other SMBH component, the outer circumbinary 
disk (iFarris et al.ll2015h . and the infalling stream. We ignore 
those complications in the fitting and will discuss the related 
effects in Appendix. 

In the standard thin accretion disk model, the emission from 
an annulus r — dr/2 —^ -f dr/2 of the disk is approximated 
by a black body radiation with an effective temperature of 


Teff(r) 


3GM.Macc 

STTCTBr'^ 



( 1 ) 


where G is the gravitation constant, ctb is the Stefan- 
Boltzmann constant, Macc is the accretion rate of the SMBH, 
and Tin is the radius of the disk’s inner edge. For the cir¬ 
cumbinary disk, rin,c = obbh/(1 + q) + Rn', for the mini¬ 
disk around the secondary SMBH, rin,c 3.5GM,^s/c^ = 
3.5rg_s, assuming a radiative ef ficiency e ~ 0.1 (correspond¬ 
ingly an SMBH spin of ~ 0.67: lYu & Lul20d8tlShankar et al.l 
120131) . Here c is the speed of light. We can then obtain the 
continuum emission for either the circumbinary disk or the 
mini-disk around the secondary SMBH as 

27rr 2/ic^cosi/A® 

^ Xi„ Df exp[hc/XkBTes{r)] -1 


where Di is the luminosity distance of Mrk231, h is the 
Planck constant, the Boltzmann constant, i the inclination 
angle, rout = lO^nu^c and /r,sf?RL(9) for the circumbinary 
disk and the mini-disk, respectively. Here we simply assume 
cosi = 0.8, the mean value for type 1 quasars. We find that a 
slight change of cos i does not affect our results significantly. 


3.3. Pseudo-continuum from Fe emission lines 

It has been shown that there are extremely strong Fell 
emission lines in the optical spectrum of Mrk 231. Thou¬ 
sands of Fe emission lines from the broad line emission re¬ 
gion blending together can form a pseudo-continuum. There¬ 
fore, the Fe emission must be included when fitting the SED 
of Mrk 231 . We use the template-fitting method introduced 
by iPhillipsI (Il977h . probably a standard practice, to treat the 
Fe emission in Mrk 231, in which the Fe spectrum of the 
narrow-line Seyfert 1 galaxy I Zw 1 is used to construct an 
Fe ll template. In the UV ban d, we adopt the Fe template 
of IVestergaard & WilkesI (120011) . In the optical band, the Fe 
template is const ructed according to the lis t of the Fe lines for 
I Zw 1 given in IVeron-Cettv et aH (120041) . We combine the 
UV and optical Fe templates to form a single template, con¬ 
volve it by a Gaussian function with a width of FWHMconv, 
and scale it by a factor of Ag to match the observations. The 
total flux of the Fe emission, /pe, is equal to Ag multiplied by 
the total flux of the template; the FWHM of the Fe lines can 
be expressed as FWHMpe = ^/FWHMj^^wi + FWHM2^„^. 
In the above algorithm, we assume that the Fe emissions in 


the UV and optical bands have the same width and there is no 
shift between the redshifts of the UV and the optical lines. We 
also assume that the ratio of the UV Fe flux to the optical Fe 
flux is fixed to be the same as that of I Zw 1. A detailed dy¬ 
namical model for the broadening and the shift of individual 
Fe emission lines is beyond the scope of our paper. 

We have also further checked the validi ty of adopting the Fe 
template of I Zw 1 by using a CLOUDY (iFerland et alJl201^ 
model, with the best fit continuum of Mrk 231 obtained below 
as the input intrinsic continuum. We And that the total flux 
of the UV Fe emission (Ipe.uv) to the total flux of the opti¬ 
cal Fe emission (/pe.opt) generated from the CLOUDY model 
3) is similar to that of IZw 1 if the ionization parameter 
is sufficiently small 0.001). A small ionization param¬ 
eter is compatible with the best-fit parameters of the BBH- 
disk accretion system, since the broad line emission region 
is bound to the primary SMBH, much more massive than the 
secondary SMBH, while the ionizing photons are provided by 
the mini-disk around the secondary SMBH. If /pe.uv/^Pe.opt 
is substantially smaller than that of IZw 1, then the FeLoBAL 
feature in Mrk 231 would be much less significant than that 
shown in Figures [3] and |2l 

4. THE OPTICAL-TO-UV SPECTRUM OE MRK 231 

The observed broadband spectrum of Mrk 231 is first 
shifted to the rest frame and then corrected for the Galactic 
extinction [i7(B — V) = 0.02213], and the results are shown 
in Figure |2] A number of fitting windows are chosen in or¬ 
der to avoid those strong emission or absorption lines, such 
as Hq,, H^, Oiii, Lya, and Nal D, which is the commonly 
adopted way to perform the continuum or SED fitting. The 
fitting windows adopted are 1140-1152, 1260-1300, 1340- 
1400, 1500-1650, 1860-2000, 3090-3145, 3740-3800, 4400- 
4700, 5080-5600,6060-6350, and 6750-6850A. A more com¬ 
plicated model may be developed by fitting each of the emis¬ 
sion lines with one or more Gaussian components and each 
of the absorption lines with Voigt profile or Gauss-Hermit ex¬ 
pansions; however, these lines will not affect the overall shape 
of the continuum, the focus of this paper, and thus we choose 
to simplify the presentation by not involving fittings to those 
emission and absorption lines. The wavelength range 3800- 
4400A is also excluded because of the contamination from 
the high-order Balmer lines. The wavelength range from 2000 
to 3000A is excluded to avoid the uncertainties in modeling 
strong Ee absorption lines as indicated by the EeLoBAL sig¬ 
natures found in the optical band. 

5. MCMC FITTING RESULTS 

We use the Markov Chain Monte Carlo (MCMC) method 
to obtain the best fit to the observational data in the above 
fitting windows and constrain the model parameters. The pa¬ 
rameters included in the continuum model are the total mass 
(M,), mass ratio (g), semimajor axis (obbh) of the BBH, the 
Eddington ratios of the outer circumbinary disk (/Edd.c) and 
the inner mini-disk (/Edd.s), the ratio of the outer boundary 
of the inner mini-disk to the mean Roche radius (/r,s), the 
scale factor Ag and convolution width FWHMconv of the Fe 
emission lines, and the extinction E-q-y due to the interstel¬ 
lar medium in Mrk 231. Here the Eddington ratio is defined 
as the ratio of the accretion rate of the circumbinary disk 
Afacc.c (or the secondary mini-disk Macc.s) to the accretion 
rate MEdd set by the Eddington limit (assuming e = 0.1), 

i.e., /Edd.c = Macc,c/MEdd{M,) = eMacc,cC^/TEdd(Vf,), 
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Fig. 2.— The optical-to-UV spectrum of Mrk231. The green curve repre¬ 
sents theUYobservat^sbytheCOSand FOS on board the HST (ai'chive 
data; IVeilleux et alJ2Q13l:ISnQith et alJ1995l) . the o bservations by the William 
Herschel Telescope (WHT) and Keck telescope ILeighlv et alJI2Q14h . The 
wavelength ranges for those different observations are marked at the bottom 
of the fig ure. The FOS data i s scaled up by a factor of 2 to match the COS 
data tsee lVeilleux et alJI2Q13h. The j ^ rey cu iwe represents the HST compos¬ 
ite spectrum of quasars <Zheng et al]ll99'^ . The continuum shows a sharp 
drop off at A ~ 4000 — 2500A and deficit of flux at A 4000 — lOOOA. 
The locations of the emission lines Ha, H/3, and Lya are also marked in the 
figure. 

/Edd.s = Mg,cc,s/M'EddiMt^s) = efface.sC^/iEdd(Vr,^s), 

LEddiM,) = 1.3 X lO46ergs-i(M,/lO®M0), and 

iEdd(M.,s) = 1.3 X lO4®ergs-i(M.,a/lO®M0). We first 
adopt those nine parameters (M„ q, obbh, /sdd.c /sdd.s, 
/r.s, FWHMconv, F'B-v) to fit the spectrum and ob¬ 
tain the best fit, and then fix the two parameters Ag and 
FWHMconv(~ 3000km s“^) at their best-fit values and ob¬ 
tain constraints on the other seven parameters (M,, q, obbh, 
/Edd,c, /Edd,s, /r.s, F^b-v)- In the fitti ng, we adopt an extinc¬ 
tion curve for SMC according to (lPeilll992h PI the Eddington 
ratios for the outer circumbinary disk and the inner mini-disk 
are assumed to be in the range from 0.1 to 1, since the stan¬ 
dard thin disk model may be invalid if the Eddington ratio is 
substantially smaller than 0.1. 

The top panel of Eigure [3 shows the best-fit model to the 
continuum emission from the BBH-disk accretion system. 
The overall shape of the observed continuum can be repro¬ 
duced by the BBH-disk accretion model. The best fit param¬ 
eters are M, = 1.5 x 1O®M0, q = 0.026, /Edd.c = 0.5, 
/Edd.s = 0.6, obbh = 590 AU= 2.9mpc, and /r,s = 0.11, 
respectivelyO According to the best-fit model, the circumbi¬ 
nary disk dominates the emission of Mrk 231 in the optical 
band; the mini-disk dominates the EUV emission; the sharp 
drop off at 4000 — 2500A is mainly due to the cut off of the 
circumbinary disk and the gap (or hole) opened by the sec¬ 
ondary SMBH, but not an extremely high extinction (see the 
top panel of Eigure O. It is also not a necessity to have dif¬ 
ferent extinctions in different bands to explain the sharp drop 
off. 

The second panel of Eigure |3 also shows the best fit to the 

^ We have checked that choosing a different extinction curve (e.g., the one 
for LMC) does not qualitatively affect our results. 

^ Note that the best fit value of /Edd,s = 0-6 means that the secondary 
SMBH accretes via a rate close to the Eddington limit. The numerical simu¬ 
lations suggest that the accretion rate onto the primary SMBH is smaller than 
that onto the secondary SMBH for a BBH-disk accre tion system with a mass 
ratio of a few to 25 percent (e.g. JFanis et al.lf2014l) . For the BBH system 
that best fits the continuum, the accretion rate of the primary disk should be 
< 0.01, which is via the advection dominated accretion flow (ADAF) mode 
iEsin et akll 19971) and radiatively extremely inefficient, and thus its emission 
can be neglected. This validates the omission of the primary disk emission in 
the fitting. 
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Fig. 3.— The optical-to-UV spectimm of Mrk231 and the model spectrum. 
From top to bottom, the first panel: the green curve represents the obser¬ 
vations as that in Figure [2] The red curve represents the best-fit model for 
the continuum emission from the BBH-disk accretion system, a combination 
of the continuum emissions from the circumbinaiy disk (cyan dashed line) 
and the mini-disk around the secondary SMBH (cyan dotted line). The sec¬ 
ond panel: the green curve represents the observational spectnam of Mrk 231. 
The red curve represents the best-fit continuum spectrum, a combination of 
the continuum emissions from the circumbinary disk, the mini-disk around 
the secondary SMBH, and the pseudo continuum by Fe emissions. The third 
panel: the green curve represents the residuals of the best fit continuum to 
the observations. The blue dotted curve represents the Fe absorption ob¬ 
tained from a model using SYNOW code, and the magenta solid curve repre¬ 
sents the FeLoBAL absorption of SDSS 1125+0029 scaled by a factor of 1.6. 
The fourth panel: the red curve represents the model spectrum by adding the 
pseudo-continuum due to numerous Fe emission lines and including the con¬ 
tribution from the FeLoBAL absorption. The black points represent the ob¬ 
servational data in those windows adopted in the continuum fitting as marked 
in this panel. 


continuum emission of Mrk 231, which includes not only the 
continuum emission from the circumbinary disk and the mini¬ 
disk around the secondary SMBH, but also the pseudo con¬ 
tinuum from numerous Fe emission lines. The Fe emission 
of Mrk 231 in the optical band is extremely high, and conse¬ 
quently the Fe emission in the UV band is likely to be high, 
too. In the observed UV s pectrum, however, there appear no 
strong Fe emission lines (IVeilleux et al.ll2M3l : iLeighlv et al.l 
l2014^ . Mrk231 is an FeLoBAL and should have very strong 
Fe absorptions in the UV band as indicated by the optical ab¬ 
sorption lines, such as Na ID, but there appear no very strong 
Fe absorption lines in the observed UV spectrum. The reason 
is that the Fe absorption dominates over the Fe emission in 
this band as detailed below. 

The third panel of Figure[3shows the residuals of the model 
spectrum at 2000-3200 A which reveals strong Fe absorp¬ 
tion lines. These Fe absorption features are quite similar to 
a 1.6-time scale up of the FeLoBAL abso rption feat u res of 
SDSS 1125H-0029 shown in Figure 4a of lHallet al.l ([200l 
the magenta curve), of which no observation at ^ 2000A is 
available. As an illustration, we further use the SYNOW code 
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Fig. 4.— Two-dimensional probability contours for different parameters of 
the BBH. The red stars represent the best-fit parameter values. Only 50, 000 
points (black) out of 50, 000, 000 simulations are plotted. The cyan curves 
represent 1-, 2-, and 3cr confidence contours from the inside to the outside, 
respectively. 


(iBranch et al.ll2002h . a fast parameterized synthetic-spectrum 
code, to generate a spectrum for the absorption features in 
the NUV band for Mrk231. We assume that only the wind 
that covers the surface of the disk facing toward distant ob¬ 
servers contributes to the absorption because of the optically 
thick disk, and this wind may launch at the inner edge of the 
circumbinary disk and/or from the mini-disk. For simplicity, 
we consider three species, i.e., Fe I, Fe II, and Mg II, and as¬ 
sume a temperature of lOOOOK, minimum velocity Vnun = 
1000 km s“^, and maximum velocity Umax = 8000 km s“^. 
The maximum velocity is on the same order of the escape ve¬ 
locity at obbh- The model absorption spectrum is shown in 
the bottom panel of Figure [3] It appears that the main ab¬ 
sorption features seen in the residuals can be roughly mod¬ 
eled, though there are still some discrepancies in details. In 
the bottom panel of Figure |3 the model spectrum (red curve) 
is a combination of the best fit to the continuum and a 1 . 6 - 
time scale up of the FeLoBAL absorption feature of SDSS 
1125H-0029 (cyan line in Fig.|3ll, which appears to match the 
observations well. 

Figures |4] and |5] show the two-dimensional probability con¬ 
tours and one-dimensional marginalized probability distribu¬ 
tions for the model parameters obtained from the MCMC 
fitting, respectively. The parameters M,, /Edd. and ubbh 
mainly determine the emission from the circumbinary disk 
and they are strongly degenerate with each other. If one of 
them could be determined by an independent method, the con¬ 
straints on the other two would be significantly improved. The 
emission from the mini-disk around the secondary SMBH is 
determined by q, /Edd,s. and /r_s, for which q and /Edd,s are 
also degenerate with each other. The constraints on these 
three parameters can be substantially improved if observa¬ 
tions at EUV (e.g., < lOOOA) are available. According to 
the marginalized one-dimensional probability distribution for 
each parameter shown in Figure |3 by the standard MCMC 
technique, we have the peak values of the model parameters 
as log(M,/M 0 ) = S.Stol, logg = -l.Stgl, log/Edd.c = 

-0.5t°:t log/Edd,s = -0.4t° J, log(aBBH/i?g) = 2.5l°;2 

(here i?g = GM./c^), logA,s = -0.5;° t E^-v = 
O.lO^Qgg, and the orbital period ~ 1.2('(g'^yr. Note these 
values for the model parameters are somewhat different from 
those from the best fit with minimum -value. 

The X-ray emission of AGNs and quasars is normally emit¬ 
ted from a corona structure above the inner disk area in 




Fig. 5.— One-dimensional probability distributions of different parameters 
of the BBH. The red dotted vertical line in each panel represents the value of 
the best fit adopted in Figs.[2and[3]to produce the model spectrum. 

the standard d isk-corona model (iHaardt & MaraschU 119931 : 
iDaiet al.l20f^ . Here in Mrk 231, the outer circumbinary disk 
is truncated at a radius substantially larger than several hun¬ 
dreds of the Schwarzschild radius of the primary SMBH, for 
which the disk is too cold and a corona probably cannot be 
established above it to emit X-rays significantly. Therefore, 
the X-ray emission can only arise from the mini-disk or its 
associated corona around the secondary SMBH in the BBH- 
disk accretion model. Since the secondary SMBH is smaller 
than the primary one by a factor of ~ l/q 38, it is almost 
guaranteed that the X-ray emission from this system is much 
weaker than that from a single SMBH-disk accretion system 
with an SMBH mass of M, p + M, g and an accretion rate 

of /Edd.cAr.Edd.c (or /Edd,sM^,Edd,c)- The best fit suggests 
that the observed X-ray luminosity of Mrk 231 in the 2-10 ke V 
band is about 1 % of the bolometric luminosity from the mini¬ 
disk, which is well consistent w ith those for normal AGNs 
and quasars (iHopkins et al.ll20^ . This solves the mystery of 
the intrinsic X-ray weakness of Mrk 231. 

6. MODEL IMPLICATION AND DISCUSSIONS 

According to the model fits, the intrinsic continuum and 
SED of Mrk 231 is significantly different from the canon¬ 
ical ones of normal quasars. The deficit of intrinsic UV 
emission will lead to substantially weaker broad line emis¬ 
sions compared with those of normal quasars. However, a 
number of broad emission lines, such as H^, H^, are evi¬ 
dent in the optical spectrum of Mrk 231. The ratio of the 
Ha flux to the H/3 flux is roughly 3, consistent with those 
of normal quasars. Here we check whether the total num¬ 
ber of ionizing photons emitted by the central source is large 
enough to balance the total number of recombinations oc¬ 
curred in the broad line region. The total number of H^ pho¬ 
tons is directly related to the number of ionizing photons as 

= IF cb(h^t) where ft/An is the cov¬ 

ering factor, a^/aB^ 1/8.5 is the number of H^ photons 
produced per hydrogen recombination, and vq = c/912 A, 
vUf) = c/4861A (lOsterbrock & Eerlandl 1200^ . Subtract¬ 
ing the best-fit continuum from the observed spectrum from 
4700A to 5800A, we obtain Wp luminosity by integrating 
the residual spectrum. Using the above equation, we find 
H/dTT ^ 0.29, which is fully consistent with the typical range 
of H for AGNs/quasars and suggests that the ionizing photon 
emission from the mini-disk around the secondary SMBH is 
sufficient to produce the optical broad emission lines, such as 
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Ha and H/3. 

The inclination angle of the disk in Mrk231 could be 
larger, e.g., cost ^ 0.5, since Mrk231 is an FeLoBAL 
(though not a normal FeLoBAL in the BBH-disk accretion 
scenario) and FeLoB ALs were suggested to h ave a larger 
i and a smaller cos* (iGoodrich & Millen 120051) . By alter¬ 
natively setting cosi = 0.5, the best fit obtained from the 
MCMC fitting suggests \og{M,/MQ) = 8.5^02’ log*? = 
log/Edd.c = -0.5l°;i log/Edd.s = -0.3l°;2, 
log(aBBH/i?g) = 2.5to;2, log/r,s = -0.5to;3, andF;B-v = 
0.071q'q 2. With this BBH configuration, the orbital period 
is 1.6^0 2 ysar and 0/47r ~ 0.5. It appears that the results 
are qualitatively consistent with those obtained by assuming 
cos i = 0.8. By relaxing the inclination angle as a free param¬ 
eter, we find that our results are not affected significantly. 

Numerical simulations suggest that the continuum emission 
from a BBH-disk accretion system may vary periodically due 
to the change of the material infalling r ate from the circumbi - 
nary disk to the inner mini-disk(s) (e.g. JHavasaki et al.l2008l) . 
though the variation for those systems with small mass ra- 
tios (less than a few percent) may not be significant (e.g., 
iFarris et ani2014ft . Several UV observations of Mrk231 have 
been obtained by Hubble Space Telescope (HST) over the 
ye ars. The UV contin uum of the HST COS observations 
by iVeilleux et al.l (1201^ is roughly co nsistent with the HST 
G160H observation bv lGallagher et"^ (l2002h . while it is per¬ 
haps a factor of two h igher t han the HST G190/270H obser¬ 
vation bv ISmith et al.l (119951). As seen from F igurel^ earlier 
observations by lUE (iHutchings & Nefflll987h show that the 
continuum (XLx) levels of Mrk231 at ISOOA ^ 2.30,3.20, 
and 3.21 x 10“^^erg cm“^ at three different epochs, and 
later on they are about 1.85 a nd 1.77 x 10~^^erg cm“^ 
measured by the HST COS (IVeilleux et al.l |20T3 ) and FOS 
(iGallagher et al.l 120021) . The probability for those observa¬ 
tions to be consistent with a constant (or no variation) is 
0.00013 according to the -statistics, which suggests that the 
UV emission ov er the past seve r al dec ades does vary, confirm¬ 
ing the claim bv iVeilleux et al.l (120131) . However, these obser¬ 
vations are not sufficient for investigating the expected (quasi- 
)periodical variations as done for the recentl y reported BBH 
candidate PG 1302-102 (iGraharn et al.ll2015h because (1) the 
total number of the observations is small; and (2) the intervals 
between some of these observations are too large compared 
with the expected orbital period. Although there are more op¬ 
tical observations of Mrk 231, the number of observations at 
each individual band is limited and not sufficient for an anal¬ 
ysis on the (quasi-)periodical variation. 

Note that t he periodical variati on of the BBH candidate 
PG 1302-102 (iGraham et al.ll20l'5h has an amplitude of only 
^ 0.14 magnitude (or ^ 14%), which suggests that the ampli¬ 
tude of the expected periodical variations of some BBH-disk 
accretion systems may be small. Future intensive monitor¬ 
ing the UV continuum emission of Mrk 231 may reveal the 
periodical variation, which would confirm the BBH hypothe¬ 
sis and be useful to further constrain the dynamical interplay 
between the BBH and its surrounding accretion flow. 

It has been shown that the polarization fraction of the 
optical-to-UV continuum of Mrk 231 depends on frequency 
and the peak of the polarization is around the wavelengths 
~ 3000A (ISmith et al.l 119951) . This dependence is proba¬ 
bly due to the scattering clo uds distributed as ymmetrically 
about the illuminating source (ISmith et al.lll995l) . In the BBH 


197B-12-22 lUE 
1980-04-13 lUE 
1980-07-26 lUE 
1996-11-21 FOS 


Fig. 6 .— UV continuum flux at 1300A (Z/ 1300 ) for Mrk231 at different 
epochs. Different points are obtained from the UV spectra obtained by lUE 
and HST at different epochs as marked in the figure, smoothed over a wave¬ 
length range of 20A around ISOOA at the rest frame. The bars associated 
with each point represent the l-cr standard deviation of the measurement. 


scenario, the blue photons around 3000A are mainly emitted 
from the inner edge of the circumbinary disk, where the disk 
may be puffed up because of the accumulation of accreting 
hot material there, which may lead to more significant scat¬ 
tering of those blue photons emitted from that region and thus 
a high polarization fraction; while those photons at the FUV 
and optical bands may experience less scattering as they are 
away from the inner edge of the circumbinary disk. This may 
explain the fr equency dependence o f th e polarization fraction , 
qualitatively: IVeilleux et al.l (120131) and lLeighlv et al.l (l2014h . 
likewise, also provided qualitative explanations on these po¬ 
larization measurements in their competing models. 

The standard accretion disk model adopted in the fitting is 
simple, and the torque due to the BBH on the outer disk is not 
co nsidered. As sh own in the Appendix, we adopt the model 
by iRaflko^ (l20 1 3h to fit the continuum, in which both the in¬ 
ternal viscosity and the external torque by the BBH on the 
circumbinary disk are considered, and we And that there are 
no significant differences in the constraints on the model pa¬ 
rameters. 

According to the best fit, the infalling rate from the cir¬ 
cumbinary disk to the inner disk(s) is smaller than the ac¬ 
cretion rate of the circumbinary disk by a factor of ^ 30. 
Currently it is not clear whether such a small infalling rate 
can be realized in BBH-disk accretion systems. For close 
BBH-disk accretion systems, a number of simulations and 
analyses suggested that the infalling rate into the gap or the 
central cavity is substantially smaller than the accretion rate 
at the outer boundary of the circumbinary disk because of 
the tidal barrier by the central BBH (e .g., the ID simula- 
tio ns by i Milosav hevic & Phinnevll2005L the 3D si mulations 
by iHavasaki et al.l 20071 or the simple arguments by iRaflkovI 
120131) . However. some recent simulations showed that the in¬ 
falling rate may not be significantly suppressed by the tidal 
barrie r (e.g., iNoble et^l2012t iMacFadven & Milosavlievia 
|2008^ and it may be even comparable to the accretion rate 
in the case with a single central MBH (IFarris et al.l 120141: 
IShi & Krolikll20L^ . Considering that all those simulations 
are quite idealized, it is also not clear whether the accretion 
onto the central BBH is really significantly suppressed or not. 
Future more realistic simulations may help to demonstrate 
whether the small infalling rate found for the BBH-disk ac¬ 
cretion system in Mrk 231 is possible. 

The spectral features of Mrk 231 were recently proposed 
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to be explained by either (1) a hybrid model with both a 
young star burst (for the UV continuu m) and an obscured 
quasar (for the optical continuum) (see iLeighlv et al.l 120141) 
or (2) an absorption model with extremely large dust red¬ 
dening at the NUV band but small reddening at the FUV 
band (see IVeilleux et al.ll201^ . In the first model, the FUV 
emission is dominated by the young star burst with an age 
of ^ 100 Myr and therefore would have no significant short 
time-scale variations, which appe ars inconsistent with the ob¬ 
servations. IVeilleux et al.l (120131) also argued that the broad 
asymmetric Lya line cannot be produced if the FUV emission 
is from an extended star burst. While in the second model, 
a small fraction of FUV emission is leaked through the ob¬ 
scuring material, and the UV variability will depend on the 
stability of the leaking holes. Future monitoring observations 
of Mrk231 at the FUV and NUV bands, and simultaneous 
multi-wavelength observations will be helpful to confirm the 
BBH-disk accretion explanation and further constrain the or¬ 
bital evolution of the BBH in the core of Mrk231. 

7. CONCLUSIONS 

In this paper, we show that various unique features in the 
optical-to-UV spectrum and the intrinsic X-ray weakness of 
Mrk 231 can all be well explained, if a pair of SMBHs exists 
in the core of Mrk 231, with the masses of the primary and 
the secondary SMBHs as ~ 1.5 x IO^Mq and 4.5 x 1O®M0, 
respectively. The existence of a BBH in Mrk 231 is compat¬ 
ible with its disturbed morphology and tidal features, which 
indicates a merger event in the past. (Note that the secondary 
SMBH is rather low in mass; however, it should be able to 
sink down to the center because the stars initial ly asso ciated 
with it enhance the dynamical friction see 120021 ) The 
semimajor axis of this BBH is ^ 590 AU, about 190 times 
of the Schwarzschild radius of the primary SMBH, and its or¬ 
bital period is just ^ 1.2 year, relatively short among the few 
known BBH candidates (IValtonen et al.l 1200^ iGraham et alJ 
120151) . which makes it an ideal system to study the dynamics 
of BBH systems. Such a BBH emits gravitational wave on 
tens of nanohertz, and the change rate of its orbital period due 
to gravitational wave radiation is about 40 seconds per orbit. 


This BBH might be a target for gravitational wave studies in 
future. 

The orbit of such a BBH system decays on a timescale 
of a few times of 10^ year due to gravitation al wave radia¬ 
tion a nd the torque of the circumbinary disk (iHaiman et alJ 
l2009h . which is not too small compared wit h the lifetime of 
quasars (a few times 10^ to 10^ year; s ee lYu & Tremainel 
l2002t IShankar et alJl20l3t I Yu & Lull2008^ . The majority of 
quasars are believed to be triggered by mer gers of galaxies 
and c onsequently involve mergers of SMBHs (IVolonteri et al.l 
120031) . and those BBH systems with mass ratio in the range 
of a few percent to 1 may lead to a notch in the optical- 
to-UV continuum emission if their semimajor axes are in 
the range of a few hundreds to a bout one thousands grav¬ 
itational radii (iRoedig et al.l 120141) . whi ch correspond to or - 
bital decay timescales of 10^-10® year (IHaiman et alJl200^ . 
Therefore, the occurrence rate of active BBH systems, with 
deficits in the optical-to-UV emission, may be rou ghly a few 
thous andths to about one percent among quasars (lYan et aP 
|2013). Our analysis of Mrk231 demonstrates the feasibility 
of finding BBH systems by searching for the deficits in the 
optical-to-UV emission among the spectra of quasars, a new 
method proposed by a number of authors (Farris et all 1201 51 : 
iGiiltekiri & MilleJl2012l : ISesana et al.ll2012 : lYan et al.ll2014l) . 


in addition to those current practices by searching for the kine- 
matic and image signatures of BBHs among AGNs/quasars 
(iBoropn & Laueii 12(3091: IValtonen et £0120081: [Grahain et alJ 


201.51 Tin et al. 2014IGaskelll 

19961 Bogdanovic et al. 

2008; 

Tsalmantza et al.l 201 It Eracleous et al.ll2012l l.lu et al. 

2013; 

Rodriguez et alJl2006llPoDoviq|2012|). 
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APPENDIX 


The structure of the circumbinary disk may be different from that of a standard thin disk, especially at th e region close t o the 
inner edge, since the torque raised by the central binary may lead to gas accumulation there. According to iRafiko^ (120131) . the 
evolution of the circumbinary disk under the actions of both the internal viscosity and the external torque can be described by a 
simple equation 


as _ 1 a a / 3 au\ 2 sa 

dt r dr \dr) \ dr) ~^ 


( 1 ) 


Here t is the evolution time, A is the external torque per unit mass of the disk due to the central BBH, E, z/, and r are the 
surface density, kinematic viscosity, and radius of the disk, respectively, and I = n{r)r‘^ is the specific angular momentum of 
gas material on a circular orbit with a radius r. At each radius, the accretion disk properties can be characterized by the viscous 
angular momentum flux 

Fj = —27rz/Er^^^ = 37rz/EUr^, (2) 

dr 

where the last equality is valid for a near Keplerian disk with U ~ and M, = M, p -f Af, g. 

Since A in equation ([T]) is significant only in a narrow annulus at the inner edge of the disk, we can assume that A = 0 outside 
of some radius ta, which is not too different from rjn. Outside ta, the disk evolves according to Equation ([T]) with A = 0. 
The so lution of this equation can be obtained by the standard approach introduced in Lynden-Bell & Pringle G974) lLvnden-BeIil 
(11974ft . if the inner boundary condition (IBC) is assumed to be torque free. However, the central BBH exerts a torque on the inner 
boundary of the circumbinary disk, which can be approximated by Fj(rin) Ri —dLBBu/dt = —TBBHfBBH/2aBBH, where 
Tbbh = / (1 + 9 )^ is the orbital angular momentum of the BBH, and ubbh = rfoBBH /dt is the inspiralling 
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speed of the BBH. According to Rafikov (2013f lRafikovl (l2013h . this IBC may be written as 


dFj 

dl 



M{rin) = xMoo, 


(3) 


where x is assumed to be a constant (< 1). By this setting, M(rin) can be substantially smaller than Moo, and disk material 
may gradually accumulate near the inner boundary if x < 1- We set the accretion rate at the outer boundary as Moo- A self¬ 
similar solution can be obtained for the evolution of the circumbinary disk once the initial separation of the BBH is given [i.e., 
aBBH(i = 0) = obbh.o]- With the solution of Fj{r, t), the effective temperature of the disk can be obtained as 


Tee{r,t) 


3 

Stt 


Fj{r,t)n 

o-Br^ 



(4) 


By integrating the multi-color black body emission over the whole circumbinary disk, we obtain the continuum flux as 

/■i-out.c 27rr 2hc^cosi/X^ 


Fx,c = 


Df exp[/ic/AA:BTeff(r)] - 1 


dr. 


(5) 


Here D\ is the luminosity distance of Mrk231, h is the Planck constant, is the Boltzmann constant, 
nn,c = aBBH(f)/(l + ?) + ?'out,c IS Set to be lO^GM./c^, and aBBH(f) is the semimajor axis of the BBH at time t. The 
evolution of obbh (t) is controlled by 

dttBBH _ ObbH OBBH ,,, 

dt few tj 

where few = [5(1 + <z)^/64g](i?g/c)(aBBH/^g)^ is the gravitation wave radiation timescale iPetersI (119641) . Rg = GM,/c^ is 
the gravitational radius, and tj = Lbbh/^Fj is the characteristic timescale for the BBH orbit shrinking due to the coupling to the 
circumbinary disk. If we adopt this model to describe the circumbinary disk and adopt the standard thin disk model to describe 
the mini-disk, we can also fit the continuum of Mrk231 by the MCMC technique. Under this approach, /Edd,s = x/Edd,c/9- 
The best fit is shown in Figure]?] and the constraint on obbh.o and t are 355i?g and 2.1 x 10^ year, respectively. Other best-fit 
parameters are M, = IO^'^Mq, /euu.c = 0.4, q = 0.02, x = 0.04 (corresponding to /Edd,s = 0.7), ubbr = 343i?g ^ 540 AU 
(corresponding to rin.c = 402i?g, /r_s = 0.12), and Eb-y = 0.14. These results are roughly consistent with those shown in 
Figure |3] which are obtained without considering the torque of the BBH on the inner boundary of the circumbinary disk. 

The torque due to the outer ci rcumbinary disk on the mini-disk around the secondary SMBH and the sh ock induced by the 
infall streams onto the mini-disk (iLodato et al.ll2009]lKocsis et al.ll20l3:lFarris et al .11201 Sc iRoedig et al.ll2014h . which are omitted 
in this study, may also introduce some errors to the fitting. Although many simulations and analysis suggest a flux deficit in the 
contin uum emission from a BBH-disk accretion system because of the gap (or hole) in the disk, including the one bviRo edig et akl 
(l2014h with consideration of the shock heating hot-spot on the inner mini-disk(s), one new simulation by iFarris et al.1 (1201.51) 
suggested that there may be no strong flux deficit in the continuum emission of a BBH-disk accretion sy stem by considering the 
emission from the infalling streams. However, only a specific case with a mass ratio of 1 is considered in IFarris et al.l (120151) and 
the stream emission therein is approximated as a thermal emission without considering inverse Compton scattering a nd radiative 
transf e r. Therefore, their res ults might not be applicable to other cases and need to be improved (see the discussion in iFarris et al.l 
1201-51) . iRoedig et al.l (I2014^ argued that the shock induced by the infalling streams on the mini-disk(s) may enhance the X-ray 
emission at ^ 100 keV, which does not affect the X-ray emission with energy significantly less than 100 keV and thus does not 
affect our estimate on the 2 — lOkeV emission. 
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Fig. 7.— The optical-to-UV spectrum of Mrk231 and the model spectrum. Legend is similar to that for Fig.|3] except that the torque at the inner edge of the 
circumbinary disk is considered in the model fitting to the continuum emission. 
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